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Destructible Bars in Disk Galaxies 
under the Dynamical Influence of a Massive Central Black Hole 
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Abstract 

The characteristics of the galactic bars that are prone to suffer a damaging impact from a massive 
central black hole are examined using flat stellar disks. We construct three disk model groups that consist 
of exponential disks with one type of velocity distribution and Kuzmin-Toomre disks with two different 
types of exact equilibrium distribution function. For each disk model group, three disks that have different 
typical Toomre's Q values are evolved to form bars through dynamical instability. Once a bar is fully 
developed, a black hole (BH), whose mass is 1% of the disk mass, is adiabatically added at the center of 
the disk. Our results show that lower-amplitude bars, that is, weaker bars are dissolved more easily by 
that BH. We have found that this destructibility is rooted in the characteristic feature that the bar formed 
spontaneously becomes shorter in length and rounder in shape with decreasing bar amplitude. Since such 
weaker bars are found to originate from colder disks in each disk model group, it follows that for a given 
form of velocity structure, the coldness of an initial disk determines whether the bar produced in that disk 
is favorable to dissolution induced by a massive central BH. In addition, the existence of bar-dissolved 
galaxies of the kind studied here is also discussed. 

Key words: black hole physics — galaxies: evolution — galaxies: kinematics and dynamics — galaxies: 
structure — methods: n-body simulations 



1. Introduction 

The existence of the supermassive black holes, whose 
masses range from ~ 10 6 M Q to ~ 10 9 Mq, that re- 
side in the centers of disk galaxies is gaining a firm 
basis observationally (Magorrian et al. 1998; Trcmainc 
et al. 2002; Marconi & Hunt 2003; Marconi et al. 2003, 
2004; Hiring & Rix 2004; Atkinson et al. 2005; Pastorini 
et al. 2007; Beifiori et al. 2009; Siopis et al. 2009). Such 
a large-mass black hole (BH) could dissolve a bar in a 
barred galaxy that embraces it, as a dynamical conse- 
quence. This kind of bar dissolution was first studied by 
Hasan and Norman (1990), who examined the stellar or- 
bits moving in combined fixed potentials of a rotating bar 
and a central BH in two-dimensional configurations. Next, 
Hasan, Pfenniger, and Norman (1993) extended that work 
to three-dimensional configurations. Both studies showed 
that central mass concentrations like massive central BHs 
can destroy a bar by converting bar-supporting orbits into 
chaotic ones. 

Although such orbital studies mentioned above are use- 
ful for understanding the mechanism of bar destruction, 
they cannot provide answers to the questions of how much 
BH mass is needed, and how long it takes, precisely to de- 
stroy a bar. As a helpful tool for resolving these issues, 
self-consistent iV-body simulations have played an impor- 
tant role in investigating the detailed evolution of a bar 
under the influence of a massive central BH. As pioneering 
work on TV-body studies, Norman, Sellwood, and Hasan 
(1996) explored the evolution of a Kuzmin-Toomre (K-T) 



disk (Kuzmin 1956; Toomre 1963) both in two- and three- 
dimensional configurations in which a bar was formed by 
the bar instability, and then a BH was added at the center 
of the disk as an external field. They showed that regard- 
less of disk geometries, the minimum BH mass required for 
bar destruction is 4 — 5% of the disk mass. Subsequently, 
Shen and Sellwood (2004) have obtained a similar BH 
mass to destroy a bar with high-quality A-body simula- 
tions using a three-dimensional K-T disk embedded in a 
halo. 

Unfortunately, K-T disks are inadequate in that real 
disk galaxies arc well-described by exponential surface- 
density distributions (Freeman 1970). Then, Hozumi 
and Hernquist (2005, hereafter HH) adopted a two- 
dimensional exponential disk and examined the effect of 
a central BH on a bar generated in that disk. Although 
the motions of stars in their simulations are restricted to 
a single plane, they have demonstrated that a BH with 
a mass of only 0.5% of the disk mass is sufficient for bar 
dissolution. This BH mass is about an order of magnitude 
smaller than that obtained by Norman et al. (1996) and 
that by Shen and Sellwood (2004). Concerning this dif- 
ference, HH argued that since exponential disks are more 
centrally concentrated than K-T disks, central BHs could 
have a more destructive impact on a bar in the former than 
in the latter, and that as a result, bars in exponential disks 
could be destructed more easily. However, Athanassoula, 
Lambert, and Dchncn (2005, hereafter ALD) have shown 
that complete bar destruction needs a BH with a mass 
of at least 5% of the disk mass even for more realistic 
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models that are composed of an exponential disk and its 
surrounding live halo. Their result suggests that the dif- 
ference in the minimum BH mass does not necessarily orig- 
inate from the difference in the mass profile of the disks. 
Consequently, the discrepancy in the minimum BH mass 
still remains unexplained, although the difference in disk 
geometry and situation between the two studies might be 
related to that difference in the BH mass: HH used a two- 
dimensional disk without a halo, while ALD adopted a 
three-dimensional disk with a live halo. 

If the BH mass required for destroying a bar is at least 
about 5% of the disk mass, bar dissolution of the kind 
stated here would be practically impossible. This is be- 
cause the corresponding BH mass amounts to approxi- 
mately 10 9 5 M Q for a typical disk galaxy, while the largest 
BH mass derived observationally from nearby spiral galax- 
ies is about 10 9 A7q (Kormendy 1988; Kormcndy et al. 
1996; Marconi & Hunt 2003). On the other hand, if that 
mass is around 0.5% of the disk mass, we can expect the 
existence of bar-dissolved galaxies in the real Universe. In 
this case, as Larson (2010) has argued, dissolving a bar 
could weaken and ultimately cease the central activities 
of disk galaxies, since a barred structure is considered ef- 
fective in fueling gas into the nucleus and this gas inflow 
could drive such activities. Therefore, also from this view- 
point, it is significant to determine the precise lower limit 
of the BH mass for bar destruction in order to reveal the 
secular evolution of the barred galaxies that harbor central 
BHs. However, it is conceivable that the BH mass needed 
to destroy a bar depends on the characteristics of the bar. 
Indeed, the bars formed by the bar instability will have 
a diversity of shapes and structures that arise from the 
difference in the kinematic structure of the disk like that 
represented by Toomre's (1964) Q profile, even though its 
mass profile is fixed. Nevertheless, in the previous jV-body 
studies noted above, the kinematic structure of the disk 
for a given mass profile was never changed systematically, 
with other parameters of the system intact. Instead, some 
parameters included in the system such as a BH mass and 
a BH softening length were surveyed systematically for a 
disk with given mass and velocity profiles. 

In this paper, we examine what type of bar is destroyed 
easily by a massive central BH, and thereby aim at un- 
derstanding what causes that difference in the minimum 
BH mass required for bar dissolution which has been de- 
scribed above. In section 2, we show disk models with 
different kinematic structures for given mass profiles in 
order to produce a wide variety of bars. In addition, a 
BH model is also provided. Our numerical method is de- 
scribed in section 3. In section 4, we present the evolution 
of the disk models, in which bars form spontaneously and 
then they are deformed by adding a BH at the center of 
the disk. The properties of the bars are disclosed as well. 
In section 5, we discuss destructible bars under the dy- 
namical influence of a massive central BH. Conclusions 
are given in section 6. 
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2. Models 

We adopt infinitcsimally thin stellar disks in order 
to avoid another damaging influence on a bar arising 
from the fire-hose instability inherent in three-dimensional 
disks (Raha et al. 1991; Debattista et al. 2004; Martinez- 
Valpuesta & Shlosman 2004). In addition, we do not in- 
clude bulge and halo components. Two different mass 
profiles for such flat disks are employed, that is, exponen- 
tial and Kuzmin-Toomre disks that are globally unstable 
to the formation of a bar. On the other hand, as an ex- 
ternal field, we model a BH that is placed at the center of 
the disk after the full growth of a bar. 

2.1. Exponential Disks 

The surface-density distributions of exponential disks, 
^ exp , are given by 

lid 

Mcxp(-R) = ^J^exp(-Rfh), (1) 

where Ma is the mass of the disk, h is the scale length, and 
R is the distance from the center of the disk. The disks are 
truncated at R = lOh. As adopted by HH, the full phase 
space is realized with the method devised by Hernquist 
(1993), in which the velocity structure is approximated 
using moments of the collisionless Boltzmann equation. 
In so doing, following the observations of disk galaxies 
(van der Kruit & Searle 1981; Lewis & Freeman 1989), 
we determine the square of the radial velocity dispersion, 
<tr 2 , so as to be proportional to the surface density such 
that 

a R 2 {R)=a Q 2 exp{-R/h), (2) 

where <jq is the radial velocity dispersion at the cen- 
ter. Although the initial models thus constructed are not 
based on exact equilibrium distribution functions, we did 
not see any violent changes in the early stages of evolution. 
This behavior indicates that these models are sufficiently 
close to exact dynamical equilibrium. 

We set up three models that have different Toomre's 
(1964) Q profiles shown in figure la. These models are 
named models El, E2, and E3, whose reference Q val- 
ues at the supposed solar radius, 2.43/i (8.5 kpc, see be- 
low), arc 1.2, 0.91, and 0.64, respectively. As explained 
by Hernquist (1993), there is difficulty in completely ful- 
filling the or 2 distribution given by equation (2) near the 
center, so that some reduction in <tr 2 is needed within a 
certain radius. As a result, the an 2 profiles of the models 
created here are not shifted vertically to one another at 
small radii. This behavior applies also to the Q profiles of 
the models, taking into consideration the definition of the 
Q parameter. Accordingly, model El has a rather differ- 
ent Q profile at small radii than models E2 and E3. Thus, 
in this disk model group, except at small radii, the func- 
tional form of the velocity structure is the same, and the 
radial velocity dispersion at each radius is different from 
model to model by a constant factor. 
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Fig. 1. Toomrc's Q parameters as a function of radius for (a) models El, E2, and E3 of the exponential disks, (b) models TK6, 
TK8, and TK10 of the Kuzmin-Toomre disks with Kalnajs's distribution functions, and (c) models TM3, TM4, and TM5 of the 
Kuzmin-Toomre disks with Miyamoto's distribution functions. Note that Q values at small radii for model El are smaller than those 
for model E2. Here, the radius is normalized by the exponential scale length h. 



2.2. Kuzmin-Toomre Disks 

The surface-density distributions of Kuzmin-Toomre 
(K-T) disks, /Ukt, are provided by 

-3/2 
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(3) 



The disks are again truncated at R = lOh. 

In contrast to the exponential disks, the K-T disks are 
known to possess exact equilibrium distribution functions 
(DFs), Fq, in analytical forms. In general, for a flat ax- 
isymmetric galaxy, the equilibrium DF of directly rotating 
stars, F + , is composed of the energy, e, and angular mo- 
mentum, j, of a star per unit mass. This prograde part of 
the equilibrium DF, F + (e,j), has been given by Kalnajs 
(1976) and by Miyamoto (1971). Here, we use both of 
the DFs. 

On the other hand, there is no definite way of prescrib- 
ing the distribution of retrograde stars. Then, we intro- 
duce them in the same manner as that adopted by Nishida 
et al. (1984). Consequently, the equilibrium DF is written 
as 



_ J (l/2)F+(e) + F+(e,j) j > 0, 



(1/W(e) 



3<0, 



(4) 



where the functions F^~(e) and F^(e,j) are derived from 
the expansion of F + (e,j) as 



F+(e,j)=F+(e)+F+(e,j). 



(5) 



The two types of DF described above involve a model 
parameter that assigns the kinematic structure to the disk. 
When denoting this parameter as tok for Kalnajs's DFs 
and as Tom for Miyamoto's DFs, we take wik = 6,8, and 
10, which are termed models TK6, TK8, and TK10, re- 
spectively, and ?7iM = 3,4, and 5, which are named models 
TM3, TM4, and TM5, respectively. We present the Q pro- 
files for models TK6, TK8, and TK10 in figure lb, and 
those for models TM3, TM4, and TM5 in figure lc. The 
models generated from Kalnajs's DFs have slightly declin- 
ing Q distributions with radius except for model TK6 that 
shows a slight increase in Q from R ~ 1.3h with radius, 



while those from Miyamoto's DFs yield steeply rising Q 
distributions with radius. At any rate, as the model pa- 
rameters, wk and Tom, increase, the Q values at all radii 
decrease for both DFs. 

In these disk models, the functional form of the velocity 
structure is completely different between these two types 
of DF. For Miyamoto's DFs, the square of the radial ve- 
locity dispersion is represented by 

$kt(R) 



<tr(R) 



2 to 



M 



(6) 



where $kt is the potential of the K-T disk. Thus, in this 
disk model group, the radial velocity dispersion at each 
radius is different from model to model by a constant fac- 
tor that depends on the model parameter Tom, so that the 
Q profiles are shifted vertically to one another, as is obvi- 
ous in figure lc. On the other hand, Kalnajs's DFs lead to 
rather complicated functional forms of the radial velocity 
dispersion. Even so, the resulting Q profiles show sequen- 
tial changes designated by the model parameter tok, as 
seen from figure lb. 

2.3. Black Hole 

A black hole (BH) model is the same as that employed 
by HH. That is, a BH is represented by a softened point- 
mass using a spline-kernel (Hcrnquist & Katz 1989), and 
its scale length is set to be O.Ol/i. 

The BH is added at time isH after a bar is fully devel- 
oped. To avoid a sudden change in the dynamical state 
of the disk caused by an addition of the BH, we increase 
the BH mass, M,(t), from to Mbh gradually with time 
t as follows: 



M.(t) 



M bh t 2 (3-2t) 0<t<1, 



Ah 



BH 



T>1, 



(7) 



where r = (t — iBH)Agrow, and t gr0 w is the time for the 
BH to attain to its preassigned mass ATbh- Thus, the 
BH is made to grow adiabatically by taking i grow to be 
sufficiently long. In our dimcnsionlcss system of units de- 



scribed in the next section, we choose t g 
set £bh to be 200. 



to be 100, and 
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According to recent observations, the Sa galaxy NGC 
4594 (Kormendy 1988; Marconi & Hunt 2003) and the SO 
galaxy NGC 3115 (Kormendy et al. 1996) are thought to 
harbor a central BH with a mass of approximately 10 9 Mq, 
which is the largest BH mass inferred in nearby spirals. 
Since this BH mass is close to 1% of the disk mass for 
a typical disk galaxy, we set Mbh = O.OIA/^- With this 
choice, we can evaluate the largest impact of a central BH 
on a bar from a realistic point of view. 

3. Method 

We use a self-consistent field (SCF) method (Hcrnquist 
& Ostriker 1992) to follow the evolution of the disks con- 
structed in the previous section and that of the systems 
including the BH. In this approach, Poisson's equation is 
solved by expanding the surface density, /z, and potential, 
<t>, of the disk in a bi-orthogonal basis set such that 



fi(R) = A nm (i)n nm (R) 



and 



n m, OR), 



(8) 



(9) 



where \i nm and $ nm are, respectively, the surface-density 
and potential basis functions, A nm are the expansion co- 
efficients, and R is the position vector. In the SCF code, 
we adopt Aoki and lye's (1978) basis set that is suitable 
for flat stellar disks. The outline of SCF simulations is ex- 
plained by Hcrnquist and Ostriker (1992) and by Hozumi 
(1997). In addition, the computational details concerning 
the simulations presented here arc described by HH. 

The disks are realized with N = 1,000,000 particles of 
equal mass. The equations of motion are integrated in 
Cartesian coordinates with a time-centered leapfrog al- 
gorithm (e.g., Press et al. 1986). We present results in 
the system of units such that G = Md = h = 1, where G 
is the gravitational constant. These units can be scaled 
to physical values appropriate for the Milky Way using 
h = 3.5 kpc and M d = 5.6 x 10 10 M o , with the result that 
unit time and velocity are 1.31 x 10 7 yr and 262 km s _1 , 
respectively. 

We set the maximum numbers of radial and azimuthal 
expansion terms in the SCF code, n max and m max , re- 
spectively, to be n max = 24 and m max = 12. In imposing 
these numbers of the expansion terms on each model, we 
save only even m-values to avoid a sizable disagreement 
between the center of the bar and the position of the BH 
caused by asymmetric features that would originate from 
odd m-valucs. Again, the computational details arc taken 
over from those stated by HH. 

Although our present simulations include a wide range 
of timescales after introducing the BH, we use a fixed 
time step that is made as small as possible. In order 
to reduce computation times, we have parallelized the 
SCF code in accordance with the prescription provided by 
Hcrnquist, Sigurdsson, and Bryan (1995) (see also Holley- 
Bockclmann et al. 2001). The time steps adopted are 
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Fig. 2. Time evolution of the bar amplitude before and after 
the addition of the black hole (solid lines) for (a) models El, 
E2, and E3 of the exponential disks, (b) models TK6, TK8, 
and TklO of the Kuzmin-Toomre disks with Kalnajs's dis- 
tribution functions, and (c) models TM3, TM4, and TM5 of 
the Kuzmin-Toomre disks with Miyamoto's distribution func- 
tions. The black hole is added at t = 200, and grows up fully 
at t = 300. The bar evolution without the black hole after 
t = 200 is also shown with dashed lines for each model. 

At; = 0.05 and At; = 0.002 before and after the addition of 
the BH, respectively. With these choices of time step, the 
relative energy error of the system after the full growth of 
the BH was, in all cases, smaller than 7.2 x 10 -6 . 

4. Results 

4-1. Evolution of Bars 

As the disks are evolved forward in time, linearly grow- 
ing two-armed spiral modes excited in the disks are de- 
formed into bars in nonlinear regimes through the bar 
instability. Wc evaluate the bar amplitude by the ab- 
solute value of the expansion coefficient, |^4 2 2(^)|, which 




E3 TK10 TM5 




;=200 f=800 r=200 f=800 t=200 f=800 

Fig. 3. Density contours in barred structures at the time just before the addition of the black hole (t = 200) and at the end of the 
simulations (t = 800) for (a) models El, E2, and E3 of the exponential disks, (b) models TK6, TK8, and TK10 of the Kuzmin-Toomrc 
disks with Kalnajs's distribution functions, and (c) models TM3, TM4, and TM5 of the Kuzmin-Toomre disks with Miyamoto's 
distribution functions. The bar major axis is aligned to the y axis. The contour levels are drawn at the 90, 80, 70, ■ ■ ■, and 30% of the 
peak amplitude on logarithmic scales. Black hole growth starts at t = 200 and ends at t = 300. All patterns rotate counterclockwise. 



corresponds to the amplitude of the fastest growing two- 
armed mode (see HH). Figure 2 shows the time evolu- 
tion of the bar amplitude for each model of the three disk 
model groups. From this figure, we find that all bars are 
fully developed at least by t = 200, and that their ampli- 
tudes are practically held constant over time to the end 
of the simulations unless the BH is added, although in 
an exact sense, some bars grow in amplitude very slightly 
with time and others show small decrements in ampli- 
tude with time, in response to a way of redistributing the 
angular momentum of a bar. We have thus chosen the 
time of adding the BH to be £bh = 200, as mentioned in 
section 2.3. We notice from figure 2, along with figure 1, 
that in each disk model group, the bar amplitude becomes 
lower as the disk is colder in the sense of typical Q values. 
Concerning the exponential disks, the bar amplitude of 
model El is larger than that of model E2, while Q values 
at small radii for model El are slightly smaller than those 
for model E2 as expressed in figure 1. This fact implies 
that the bar amplitude would not be sensitive to the inner 
velocity structures but to the relatively outer ones. 

Once the BH is added at t = t-Q-R , the bar amplitudes 
start to decrease with time, as demonstrated in figure 2. 
After the full growth of the BH, the bar amplitudes con- 
tinue to decay almost exponentially with time till the end 
of the simulations. Consequently, the bars become very 
round in the final stages as found from a comparison be- 
tween the density contours at t = <bh and those at the 
end of the simulations, which are illustrated in figure 3. 



To quantify to what degree the bars are made round, we 
evaluate the axis ratios of the bars along radius at these 
two times by calculating the principal moment of iner- 
tia tensor for particles included in a specified radius, the 
value of which is regarded as the axis ratio at that radius. 
Resulting axis-ratio profiles are presented in figure 4. We 
can see from this figure, together with figure 3, that in 
each disk model group, the bar is more easily deformed 
into a rounder shape as the smallest axis ratio along ra- 
dius is larger, that is, the bar is rounder when the BH is 
added, or as the bar half-length, estimated here to be a 
distance from the center to the radius of the smallest axis 
ratio, is shorter. 

Since the final axis ratios of the bars exceed about 0.9 
at all radii except for the bar of model TM3 (the final axis 
ratio is 0.86), irrespective of the mass and velocity profiles, 
it follows that a BH with a mass of 1% of the disk mass can 
destroy almost completely, or in some cases appreciably, 
the bars formed by the bar instability. Furthermore, when 
we take into account the Q profiles shown in figure 1, the 
runs of the final axis ratios in figure 4 indicate that in each 
disk model group, the bars originating from the colder 
disks are destroyed more easily by the BH. 

4-2. Characteristics of Bars 

We describe the characteristics of the bars before the 
BH acts on them to unravel what type of bar is more 
destructible under the dynamical influence of a massive 
central BH. 
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Fig. 4. Axis ratio of bars as a function of radius for (a) 
the exponential disks, (b) the Kuzmin-Toomrc disks with 
Kalnajs's distribution functions, and (c) the Kuzmin-Toomre 
disks with Miyamoto's distribution functions. Dashed lines 
show the axis ratios of the bars just before the black hole is 
added (t = 200), while solid lines represent those at the end 
of the simulations (t = 800). 

Figure 5 presents the surface-density profiles of the bars 
along the major and minor axes at t = teH just before the 
BH is added. Except for models El and TK6, there is 
no noticeable change in the density slope on logarithmic 
scales from the center towards the bar ends along the ma- 
jor axis. On the other hand, models El and TK6 show a 
bend in the surface-density distribution along the major 
axis at around an intermediate position between the cen- 
ter and the bar end, although the bend for model TK6 
is rather gentle as compared to that for model El. The 
surface-density profiles along the minor axis for these two 
models are close to those of the rest of the models. 

To examine the structures of the bars in detail, we 
analyze the surface-density profiles by calculating their 



Fourier components along radius, I m (R), where to = 2,4,6, 
and 8. These Fourier components arc normalized by 
the to = component, Iq, of the disk being analyzed. 
The resulting relative Fourier amplitudes of each model 
at t = £bh are displayed in figure 6. Here, we calcu- 
late the Fourier components in the same manner as that 
done by Ohta, Hamabe, and Wakamatsu (1990) and by 
Athanassoula and Misiriotis (2002). This figure indicates 
that for the hottest models in each disk model group from 
a viewpoint of typical Q values, that is, for models El, 
TK6, and TM3, the amplitudes for the to = 4 and to = 6 
components are large in magnitude as well as for the to = 2 
component, while for the rest of the models they are small, 
and, in particular, the amplitudes for the m = 6 compo- 
nent are negligibly small. In addition to this property 
of the hottest models, taking into account the Q profiles 
presented in figure 1, together with figure 5, we can see 
that a bend in the surface-density distribution of the bar 
along the major axis is produced for the models that have 
relatively high Q values near the center as for models El 
and TK6. 

To further characterize the bars, we measure the size of 
the bar region at t = £bh, following Ohta et al. (1990), who 
defined the bar region to be the region where the value of 
(la + h + h + h)/{h - h + h - h) exceeds 2.0. In figure 
7, the measured size of the bar region (bar half-size) is 
plotted against the bar amplitude that is normalized by 
the amplitude of the ring mode (m = 0), A 2 2|/|^oo I, at 
t = ^bh hi order to be compared among the disk models 
with the different mass and velocity profiles. This figure 
reveals that the bar size is approximately proportional to 
the normalized bar amplitude. 

As another ingredient for featuring the bars, the small- 
est axis ratio along radius at t = tBH, as a measure of the 
degree of elongation, is plotted against the normalized bar 
amplitude in figure 8. We can find from this figure that 
on the whole, the roundness of the bar decreases mono- 
tonically with increasing normalized bar amplitude. 

5. Discussion 

We have obtained a wide variety of bars from the ex- 
ponential and K-T disks to which various velocity distri- 
butions are assigned. The amplitude of these bars con- 
tinues to decay exponentially with time till the end of 
the simulations after the full growth of the BH that is 
added at the center of the disk. When we pay attention 
to the three models of each disk model group, the final 
axis-ratio profiles in figure 4 show that the BH can make 
the bar rounder as the bar amplitude becomes lower. To 
be sure, it may be no wonder that lower-amplitude bars 
are destructed more easily by a given mass BH. The im- 
portant thing is that this fragility of low-amplitude bars 
originates from the bar properties exposed in figures 7 and 
8 that the normalized bar amplitude is related to the size 
and shape of a bar. Therefore, the precise explanation is 
that since lower-amplitude bars are shorter in length and 
rounder in shape, a given mass BH can provide a more 
destructive impact on shorter bars that lie more deeply 
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Fig. 5. Surface-density profiles along the major (solid lines) and minor (dashed lines) axes at the time just before the black hole is 
added (t = 200) for (a) models El, E2, and E3 of the exponential disks, (b) models TK6, TK8, and TK10 of the Kuzmin-Toomre 
disks with Kalnajs's distribution functions, and (c) models TM3, TM4, and TM5 of the Kuzmin-Toomre disks with Miyamoto's 
distribution functions. 




Fig. 6. Fourier amplitudes of the surface-density distribution, I m , for the m = 2 (red lines), 4 (blue lines), 6 (green lines), and 8 
(black lines) components normalized by the m = component at the time just before the black hole is added (t = 200) as a function 
of radius for (a) models El, E2, and E3 of the exponential disks, (b) models TK6, TK8, and TK10 of the Kuzmin-Toomre disks 
with Kalnajs's distribution functions, and (c) models TM3, TM4, and TM5 of Kuzmin-Toomre disks with Miyamoto's distribution 
functions. 
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Fig. 7. Bar half-size as a function of bar amplitude normal- 
ized by the amplitude of the ring mode excited in each model 
disk. The bar half-size and the amplitudes of the bar and 
ring modes are calculated at the time just before the addition 
of the black (f = 200). The definition of the bar half-size is 
described in the text. 
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Fig. 8. Axis ratio of bars as a function of bar amplitude nor- 
malized by the amplitude of the ring mode excited in each 
model disk. The axis ratio means the smallest axis ratio along 
radius for each model. The axis ratio and the amplitudes of 
the bar and ring modes are calculated at the time just before 
the addition of the black hole (f = 200). 

inside the BH's "sphere of influence", and can increase 
the roundness more easily for rounder bars over a given 
period. 

We should mention the bar properties extracted in 
the present study from an observational point of view. 
Elmegreen et al. (2007) have found from if s -band images 
of barred galaxies that longer bars have higher bar am- 
plitudes. On the other hand, analyzing near-infrared im- 
ages of nearby barred galaxies, Menendez-Delmestre et al. 
(2007) have indicated a weak tendency that the bar el- 
lipticity is larger for longer bars. Here, the bar elliptic- 
ity is defined as 1 — b/a, where b/a is the axis ratio of 
a bar. Also, very recently, Hoyle et al. (2011) have pre- 
sented a similar correlation between the bar ellipticity and 
bar length from an analysis of bar length measurements 
in local disk galaxies. Then, if we are allowed to com- 
bine these two kinds of observed results, we are led to 
the conclusion that the bar ellipticity increases with in- 
creasing bar amplitude. In fact, the bar ellipticity is often 



used as a measure of the bar strength: larger bar cllip- 
ticitics mean higher bar amplitudes (Chapclon, Contini, 
& Davoust 1999; Laine et al. 2002; Laurikaincn, Salo, 
& Rautiaincn 2002). Therefore, the properties that we 
have uncovered for the simulated bars are consistent with 
those of real barred galaxies. Although the correlation of 
the bar length and ellipticity with the bar amplitude has 
thus far been explained from a viewpoint of secular evo- 
lution (Elmegreen ct al. 2007; Menendez-Delmestre et al. 
2007; Hoyle et al. 2011), we will demonstrate in a future 
paper that these bar properties are closely related to the 
properties of the fastest growing, globally unstable two- 
armed modes excited in stellar disks, which finally lead to 
bars in nonlinear regimes. 

We have also found that the lower-amplitude bars are 
generated from the colder disks in each disk model group, 
although the physical origin of this finding is unclear at 
the moment (we will mention it also in that future pa- 
per). Consequently, once the functional form of the ve- 
locity distribution is assigned, the degree of the coldness 
of a disk determines the fragility of a bar under the dy- 
namical influence of a massive central BH at least for 
two-dimensional disks. If this property is applicable to 
bars formed in three-dimensional disks, the discrepancy 
that HH's minimum BH mass necessary for bar dissolu- 
tion in the exponential disk is about an order of mag- 
nitude smaller than ALD's one might have arisen partly 
from the difference in the coldness, more precisely, the 
Q profile of the disk. However, it is uncertain whether 
the relation between the coldness of a disk and the ampli- 
tude of a bar holds for bars produced in three-dimensional 
disks, because the fire-hose instability intrinsic in such 
disks (Raha et al. 1991; Dcbattista et al. 2004; Martinez- 
Valpuesta & Shlosman 2004) and bar-halo interactions 
(Athanassoula 2002; Athanassoula & Misiriotis 2002), 
both of which are not included in the present study, might 
distort those properties of the bars which have been ob- 
tained in flat disks. Thus, we need to investigate the va- 
lidity of this relation using realistic disk galaxy models in 
three-dimensional configurations. 

Note that the relation between the coldness of a disk 
and the fragility of a bar mentioned above is obtained by 
ignoring the effects of gas components on the bar insta- 
bility, because our simulations are based on pure stellar 
dynamics. In general, if gas components are added to a 
stellar disk, the system will become colder on account of 
their low velocity dispersions. Then, we could expect a 
lower-amplitude bar in such a system that would suffer a 
more violent bar instability. In a case where a large gas 
mass fraction is included in a disk, Shlosman and Noguchi 
(1993) demonstrated that the bar instability in the disk 
is damped heavily. Even in this situation, it can be inter- 
preted apparently that an extremely weak bar or bar-like 
structure is formed in an exceedingly cold disk owing to a 
large amount of gas. Thus, in principle, our finding that 
colder disks lead to lower-amplitude bars might hold even 
for the disks that contain gas components, regardless of 
their amount. 

Figures 5 and 6 have indicated that robust bars against 
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the dynamical influence of a massive central BH have 
relatively larger Fourier amplitudes for the m = 4 and 
m = 6 components (models El, TK6, and TM3), and 
that for the models whose Q values are somewhat large 
near the center (models El and TK6), the bars repre- 
sented by these Fourier components exhibit a bend in the 
surface-density profile along the major axis. These two 
features appear similar to those of the MH-type models, 
employed by Athanassoula (2002) and by Athanassoula 
and Misiriotis (2002), whose halos are characterized by 
small core radii. On the other hand, the surface-density 
profiles of the models other than models El and TK6 have 
no noticeable change in density slope along the bar ma- 
jor axis. This feature corresponds to that obtained from 
the MD-type models, also used by Athanassoula (2002) 
and by Athanassoula and Misiriotis (2002), whose halos 
are depicted by large core radii. ALD have reported that 
for a given mass central BH, bars formed in the MH-type 
models are less destructible than those in the MD-type 
models. Interestingly, ALD mentioned that the disks in 
the MH-type models have Q = 1.4, while those in the MD- 
type models have Q = 1. Thus, this relation between the 
degree of bar destruction and the typical Q value is in 
concord with our claim that bars formed in colder disks 
are more destructible. It is true that we cannot compare 
directly between their results and ours, since the disks in 
their MH- and MD-type models are three-dimensional and 
their halos are self-gravitating, while our models are two- 
dimensional and have no halos. Nevertheless, our models 
seem to reflect the characteristics possessed by the bars in 
their MH- and MD-type models. This finding may imply 
that the velocity structure of the disk affects more sub- 
stantially the characteristics of the bars formed by the bar 
instability than the difference in the central concentration 
of self-gravitating halos. Therefore, we should examine 
the bar properties for both MH- and MD-type models that 
have a wide variety of velocity structures in the disks, and 
evaluate the minimum BH mass needed to dissolve the 
bars formed spontaneously. This line of investigation will 
enable us to confirm whether the discrepancy in the min- 
imum BH mass between HH's study and the others is due 
to the difference in the velocity structure of the disk. 

HH argued that the destructive effect of a central BH 
on a bar could be considered stronger for more centrally 
concentrated disks. Although Norman et al. (1996) used 
a K-T disk and ALD adopted an exponential disk, both 
studies resulted in a similar BH mass that is about 5% of 
the disk mass as the minimum value required for bar disso- 
lution, in spite of the fact that exponential disks arc more 
centrally concentrated than K-T disks. This result might 
suggest that the velocity structure is more significant to 
bar dissolution than the degree of mass concentration of 
the disk, because our simulations have shown that the vul- 
nerability of a bar to a given mass central BH depends on 
the velocity structure of the disk. Furthermore, we have 
found from figures 7 and 8 that the size and roundness 
of a bar arc determined by the normalized bar amplitude, 
irrespective of the mass profile of the disk, and that the 
strength of this amplitude is connected to the destructibil- 



ity of a bar. However, we cannot answer exactly whether 
these bar properties depend on the degree of mass con- 
centration of the disk, since it is prohibitively difficult to 
construct such models that have the same velocity struc- 
ture for different mass profiles. Therefore, we cannot con- 
clude definitely that the bars generated in more centrally 
concentrated disks are destroyed more easily by a given 
mass central BH. 

We have found that the bar amplitude continues to de- 
crease almost exponentially with time for 500 time units 
after the full growth of the BH. This simulated period 
corresponds to about 7 Gyr, if the time unit is repre- 
sented by the physical value appropriate for the Milky 
Way. We have also found that the bars formed in the 
exponential disks have been practically destroyed within 
that period by the BH whose mass is 1% of the disk 
mass. This BH mass is reduced to be approximately 
10 9 M Q for a typical disk galaxy, and so, it is compara- 
ble to the largest mass for central BHs in nearby spi- 
rals (Kormendy 1988; Kormendy et al. 1996; Marconi 
& Hunt 2003). Therefore, our simulations suggest that 
bar dissolution induced by a massive central BH could oc- 
cur in the real Universe, if galaxy disks are relatively cold 
at birth. However, the frequency of such bar dissolution 
would be rare, since the minimum BH mass examined by 
HH (about 10 8 5 M Q ) as well as the BH mass in this work 
are close to the upper bound of the observationally de- 
rived BH masses. Consequently, the bar fraction would 
be declined to a negligible degree, if any, as the Universe 
evolves, by the bar dissolution mechanism studied here. In 
fact, Elmegreen, Elmegreen, and Hirst (2004) and Jogee 
et al. (2004) found that the observed bar fraction is nearly 
constant, and recently, Sheth et al. (2008) have shown that 
it even increases with time over the last 7 Gyr (see also 
Abraham et al. 1999; van den Bergh et al. 2002). 

6. Conclusions 

We have examined the properties of the bars formed by 
the bar instability in the exponential disks with approxi- 
mate equilibrium velocity distributions and the Kuzmin- 
Toomrc disks with two different types of exact equilibrium 
DF in two-dimensional configurations, and the damaging 
impact of a massive central BH on these bars. Then, we 
have found that lower-amplitude bars, that is, weaker bars 
are destroyed more easily by a given mass central BH, 
because the bar formed by the bar instability is shorter 
in length and rounder in shape as the bar amplitude de- 
creases, thereby being more susceptible to the dynamical 
influence of that BH. 

In addition, we have also found that lower-amplitude 
bars are generated from colder disks for a given functional 
form of the velocity distribution. This finding suggests 
that destructible bars under the dynamical influence of a 
given mass central BH are determined by the kinematic 
properties of disks, like those specified by typical Q values, 
at the formation epoch, once the functional form of the 
velocity distribution is assigned. Thus, it is likely that 
the discrepancy in the lowest BH mass required for bar 
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dissolution found in the previous simulations might have 
resulted from the difference in the velocity structure of 
the disk. 

When scaling the dimcnsionlcss units to appropriate 
physical values, we have shown that the bar amplitude 
continues to decrease almost exponentially with time for 
about 7 Gyr after the full growth of a BH, and that the 
bars formed in the relatively cold exponential disks have 
resulted in practical dissolution for a BH with a mass of 
approximately 1O 9 M0, a likely largest BH mass derived 
observationally from nearby spirals. We thus infer that 
bar dissolution induced by a massive central BH could oc- 
cur in the real Universe, if galaxy disks are formed in a 
rather cold state. 

The author is indebted to Dr. K. Nitadori for optimiz- 
ing the two-dimensional SCF code for parallclization. He 
would also like to thank the anonymous referee for the 
valuable comments that have called the author's atten- 
tion to the connection of the present work to observations, 
which have helped improve the manuscript. Numerical 
computations were carried out on the PC cluster and 
Cray XT4 at the Center for Computational Astrophysics 
(CfCA), the National Astronomical Observatory of Japan. 
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